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I 

! Abstract 

The  idealized  problems  of  (a)  an  infinitesimal  Hertzian 
dipole  in  and  over  a perfect  dielectric  coating  a perfect 
conductor  and  (b)  as  Abraham  dipole  lying  on  the  conductor 
are  treated.  Unintegrated  forms  of  the  Hertz  potentials  are 
obtained  for  both  electric  and  magnetic  dipoles. 

Integrated  far-zone  forms  of  the  potentials  and  fields 
are  obtained  for  electric  dipoles  by  means  of  asymptotic 
integrations.  Par-zone  radiation  patterns  are  given  in 
order  to  indicate  the  distortions  of  the  fields  and  the 
magnitudes  of  the  residue  waves  caused  by  the  dielectric 
coatings. 

It  is  proved  that  the  power  radiated  by  the  dipole  may 
be  divided  into  two  independent  quantities -the  power  fed  to 
radiation  type  and  that  fed  to  surface  or  guided  type  fields. 
For  certain  cases  numerical  results  are  given  for  the  total 
power  radiated  and  the  relative  powers  fed  to  the  two  types 
of  fields. 

Formulas  are  derived  and  illustrated  with  numerical 
examples  of  the  radiation  resistances  of  the  dipoles  and 
the  attenuation  constants  of  the  surface  modes  due  to  finite 
conductivity  of  the  ground  plane. 


I. 

Introduction 

A study  of  the  coated  conductor  problem  was  undertaken  in 
an  effort  to  explain  certain  characteristics  of  an  antenna 
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field  pattern  obtained  from  measurements  taken  over  an 
aluminum  ground  screen.  The  effect  of  the  finite  conductivity 
of  the  aluminum  did  not  yield  an  adequate  explanation  of  the 
observed  characteristics  near  the  ground  screen  and  subsequent 
measurements  showed  that  the  phenomenon  could  be  emphasized  by 
the  addition  of  a dielectric  layer  to  the  ground  plane.  Since 
a dielectric  coated  metallic  surface  can  support  a guided  wave 
for  any  thickness  of  coating,  it  was  deduced  that  the  thin 
coating  of  aluminum  oxide  present  on  the  ground  screen  was 
responsible  for  the  behavior. 

The  idealized  problem  involves  a uniform  layer  of  perfect 
dielectric  covering  an  infinite,  perfectly-conducting  plane$  a 
Hertzian  dipole  is  in  one  of  three  positions  § Case  A,  above 
the  dielectric  $ Case  B,  in  the  dielectric $ or  Case  C,  an  Abraham 
dipole  on  the  plane,  Fig.  I.  Unintegrated  expressions  are 
derived  for  the  Hertz  potentials  in  and  above  the  dielectric 
due  to  both  vertical  electric  and  magnetic  dipoles.  For  the 
special  cases  of  vertical  electric  dipoles  in  or  near  fairly 
thin  dielectric  coatings  the  integrated  forms  of  the  far-zone 
fields  are  obtained.  Marcuvitz^  has  derived  a similar  ex- 
pression for  the  Hertz  potential  above  the  dielectric  due  to 
a dipole  above  but  near  a thin  dielectric  coating  by  a some- 
what different  method.  The  mathematical  details  of  this 

problem  are  similar  to  those  of  the  vertical  dipole  over  a 

2 q 

conducting  earth. 

The  resulting  fields,  comprising  a radiated-type  wave 
(or  compensating  wave)  and  one  or  more  guided-type  waves, 
have  been  combined  at  typical  distance  for  a polystyrene  coating 
to  yield  field  patterns.  Several  examples  of  these  are  given 
illustrating  the  distortion  of  the  patterns  due  to  dielectric 
layers.  To  further  illustrate  this  distortion,  expressions 
are  derived  for  the  total  power  output  of  the  dipole,  the  ratio 
of  the  powers  in  the  guided  and  : adiated  waves,  and  the  radiation 
resistances  of  the  dipoles.  *here  possible  in  the  illustrative 
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curves  of  these  expressions,  comparison  curves  for  metal 
surfaces  without  dielectric  coatings  are  given*  In  addition, 
the  expression  for  the  attenuation  constants  of  the  guided 
waves  due  to  finite  but  large  conductivities  of  the  conducting 
plane  is  derived.  No  attempt  is  made  at  this  time  to  -compute 
the  attenuation  due  to  conductivity  of  the  dielectric -►  A good 
approximation  to  this  attenuation  is  given  by  Attwood.4- 


II 

Formulation 

Let  an  infinitesimal  electric  Hertzian  dipole,  a magnetic 

dipole,  or  an  Abraham  dipole  be  oriented  along  and  parallel  to 

the  a-axis  of  the  cylindrical  coordinate  system  r,  0,  z,  Fig.  I. 

This  system  is  independent  of  0.  Region  I,  d — z ^oo , with 

constants  e,  u,  or  = o,  is  the  space  above  the  dielectric $ region 

2,  with  constants  n2e,  jx,  o*  - o,  o s£  z ^ d,  is  the  dielectric 

region^  and  region  3,  z o is  a perfect  conductor  with  o- = cq 

The  dielectric  constant  of  medium  2 relative  to  that  of  medium 
2 

1 is  n , where  n Is  the  index  of  refraction.  $ 

As  Is  well  known  the  electromagnetic  field  excited  by  electric 
or  magnetic  sources  oriented  in  the  z - direction  of  a cylindrical 
coordinate  system  may  be  derived  from  the  s - components  of  the 

- i 

electric  or  magnetic  Hertz  vectors.  Hence,  the  problem  will  be 
formulated  in  terms  of  these  components. 

The  problem  is  subdivided  into  five  divisions  according  to 
the  positions  of  the  dipoles  and  as  to  whether  the  electric  or 
magnetic  case  is  treated.  The  harmonic  time  dependence  e“iu>t  is 
assumed  throughout  and  is  cancelled  in  all  expressions. 


To  convert  to  dependence  e 


ja>t 


substitute  3 for  -i. 
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Vertical  Electric  Dipoles 


Case  A A Hertzian  dipole  above  the  dielectric  at 
d ^ z * < 00  , r 1 =0. 

2 2 p1(z»)B(z-z*)S(p) 

1.  V ni  + 3 ni  2ner Region  1 

(A)  (A) 


2*  V2n2  + n232fT2  = 0 

(A)  (A) 


Region  2 


Case  B A Hertzian  dipole  in  the  dielectric  at 


0<z'<  d9  r'  = 0. 


3 o V 2tt1  + P2"!  = 0 

(B)  (B) 


Region  1 


2 2 2 - P2^Z*  )&(z-z')&(r) 

4«  v + n P n?  p 

(B)  2ttn_er  Region  2 

Case  C An  Abraham  dipole  on  the  ground  plane  at 


z'  = 0,  r*  = 0. 


5«  V 2tti  + P2ni  = 0 

(C)  (C) 


60  V u,  * n 3 


P2(o)5(z)B(r) 


(C)  (C) 


2 

4-rm  er 


Region  2 


The  appropriate  boundary  conditions  ares 

dTT,  dTT2 

?•  a)  (r9d,z«)  = -rf  (r,d,z«) 


2 

eTT1(r,d9z,9)  * n,  engCi^d^*) 
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dn0 

c)  ZT  (r»°iz,)  = 0 

e *6 

d)  The  Sommerfeld  condition  of  radiation,  ’ 

where? 


TT^Cr^z1)  is  the  z-component  of  the  electric  Hertz 
vector  in  region  1» 

tT2(r,z,z')  is  the  z-component  of  the  electric  Hertz 
vector  in  region  2. 


These  are  defined  by: 


8. 


ardz 


d* 


© 


The  upper  quantities  in  the  brackets  are  valid  in  region  1 
and  the  lower  quantities  in  region  2. 


3 *=  a) 

o 

n * relative  dielectric  constant  of  region  2 

with  respect  to  region  1. 


P^z') 

P2(z«) 

4 

P2(°)/2 


is  a vertical  Hertzian  dipole  at  z* 


is  an  Abraham  dipole  at  z*  = 0 


d < z*  < oo 
0 < z*  < d 


The  Laplacion  operator 
in  the  cylindrical  co- 
ordinate system  with  no 
0 dependance. 


is  the  Dirac  delta-function  for  a 
at  z *=  z‘,  r = 0. 


S(z-z 1 )S(r) 
2nr 


source 
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Vertlcal  Magnetic  Dipoles 

Case  A A magnetic  dipole  above  the  dielectric  at 

11. 

12. 

Case  B 

13. 

14. 

One  should  note  that  the  problem  of  a magnetic  dipole  placed 
at  z*  =0  need  not  be  formulated  since  according  to  the  theory 
of  images  it  yields  fields  which  are  identically  zero. 

The  appropriate  boundary  conditions  are 


a) 

nml(r*d*z,)  = 

b) 

-§S"ml(r*d’z'> 

= h V>(r’d’z) 

subject  to  = p,  2 = P 

c) 

t^2(rf°,z*)  = 

0 

d) 

The  radiation 

condition. 

d < z ' <_cx>  ? r 1 =0 

P jim1(z,)6(z-z,)5(r) 

V Vl  + P "ml  = - — 2ttt Region  1 

(A)  (A) 

V ^nm2  + n2p2nm2  = 0 Region  2 

(A)  (A) 

A magnetic  dipole  in  the  dielectric  at 
0<  z'<d,  r*  = 0 

v2nml  + p2nml  = 0 Region  1 

(B)  (B) 

p p p umpCz’J&Cz-z'^Cr) 

V trm2  + n 3 rrm 2 = - 55 Region  2 

(B)  (B) 
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IWr.z.z')  1 


Where: 

‘'ml 

nm2(r»z»z,)  J 

These  are  defined  by 


is  the  z- component  of  the  magnetic 

Hertz  vector  in  Jre6l°n  1 

(region  2 


16, 


17  < 


E 


2m 


A d 

-1*  h 


n 


ml 

Tm2 


>l’2  (rf+p2{n|j 

Bm2  H 


n 


ml 


n. 


m2 


+ r r 

r dr3z 


ml 

Tm2 


Where  the  subscript  m indicates  quantities  in  the  magnetic 
cases  and 


n^Cz'  ) 
m2(z'  ) 


f is  a vertical  magnetic  dipole  at  zV 


dSz’^co 
o<z' £ d 


Since  the  problem  is  formulated  in  cylindrical  coordinates 
with  0 - symmetry  it  is-  convenient  to  apply  the  well-known 

7 

Fourier  - Bessel  or  Hatikel  Transform  pair  Tor  0 symmetry  to 
equations  1-6  and  11-14. 

18.  TrCX^z^z1  ) = rTT(r,z,z' )JQ(\r)dr 


19. 


Tr(r,Z,Z’  ) 


\n(X.,z,z’  )Jo(Xr)dX 


After  applying  equation  18  to  the  equations  1-6  and  11-14, 
performing  an  integration  by  parts,  and  applying  the  radiation 
condition*  to  the  left  sides  of  these  equations,  the  following 

?The  quantity  rJo(\r)  -nrAJ^CXr)  J appears  after  the  inte 
gration  by  parts.  This  disappears  identically  at  r=o  and  at  r 
the  radiation  condition  5,6  requires  its  disappearance. 
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results  are  obtained  for  the  left  sides  of  these  equations? 

for  equations  1,  3,  5,  11 9 13 

IT 5 'i2)  °*  “j, 


(■ 


for  equations  2,  4,  6,  12,  and  14 

(-4  - »2)  S,  „ 

az2  2 ®2 


20. 


21, 


to  2 

and  m are  defined? 


{2 


m2  = \2  - n2p2 


After  noting  that  the  application  of  equation  18  to  the 
right  sides  of  equations  1,  4,  6,  11  and  14  yields  the  value 
l/2n,  the  following  transformed  equations  are  obtained? 


22. 


6 


-62 


n. 


1 

(A) 


rr. 


1 

(B) 


nl 
(C) 

nml 

(A) 

"ml 

(B) 


P1(z,)5(z-z’ ) 
2ne 


)8(z-zl  ) 


2ir 

0 
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These  are  the  transformed  forms  -of  equations  1,  3,  5,  11,  13 
and  2,  4,  6,  12,  and  14  respectively.  These  satisfy  boundary 
conditions  7 a-d  and  15  a-d  as  can  be  seen  from  equations  18  and 
19  assuming  that  the  order  of  integration  and  differentiation, 
with  respect  to  z,  may  be  interchanged. 

According  to  the  physical  considerations  of  these  problems 
it  is  convenient  to  substitute  the  following  forms  in  equations 
22  and  23.  If  these  satisfy  the  equations  and  boundary  conditions 
they  are  the  unique  solutions  of  the  problems  except  for  an  arbi- 

O 

trary  constant  in  time.  The  subscripts  indicate  the  equation 
to  which  each  trial  solution  belongs. 


t -t\  z-zA  -f(z+z»-2d)+2md 

= e + e 


TTmia^z*; 

1*1- 

VMjiz*)  J 


+ C(\,z»)e 


rfCz+z'  -2d) 


Where  the  first  term  is  the  incident  wave  or  the  wave  radiated 
by  an  isolated  source,  the  second  term  is  a wave  reflected  from 
the  metal  surface,  and  the  third  term  represents  a wave  reflected 
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from  the  dielectric  surface 

■\ 

25.  eTT2(A.,z,z' ) 


LAI 


P^Z^) 


LA) _ 

|Uk(z*> 


f * D(\,jb.* ) e 


-[£  d-z*  + m(d-z)  3 


J 


+E(X.,z»  )e 


-[£  |d-z*  +m(d+z)] 


Where  the  first  term  is  a wave  transmitted  into  the  dielectric 
and  the  second  term  is  a wave  reflected  from  the  metal  surface. 

*\ 

2 — 

26.  n en2(X,z,z*) 


m. 


2n28TT2(X.,Z,0) 

. - (cj 


Lx,  -miz-*‘i 

' 4nm 


Po(°) 


iBi 


\m2(zt  j 


, -m(2d-z-z') 

+ F (\z  )e 


-m(z+z* ) 


+ G (*.,z*)  e 


z*  = 0 for  case  C. 


Where  the  first  term  is  the  incident  wave,  the  second 
term  is  a wave  reflected  from  the  boundary  z = d,  and  the 
third  term  is  a wave  reflected  from  the  metal  boundary. 


I 
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27  • r^en^^ZjZ.' ) 

-p*i'£n — 
2 

2n  en^(X.,z,o) 


— 

Vl(Vz,z') 


^(z7!  J 

where  this  term  represents  a wave  radiated  from  the 
dielectric  surface. 

It  mav  be  observed  that  the  left  sides  of  equations  (24-27) 
are  the  Green's  functions  for  equations  (22-23).  The  proper 
choice  of  the  coefficients  of  the  first  terms  on  the  right 
sides  of  equations  (24)  and  (26)  gaurantees  that  these  equations 
satisfy  the  conditions  for  the  Greens  functions  at  z - z'. 

Substituting  equations  equations  (24-27)  in  their  proper 
boundary  conditions  7a -d  or  15a-d,  it  may  be  verified  that  the 
trial  solutions  24-27  indeed  do  satisfy  the  boundary  conditions. 
Evaluating  the  constants  B(\,’z')  through  H(A.,z')  yields  the 
transformed  Hertz  potentials.  Applying  the  Fourier-Bessel 
integral,  19,  rearranging  the  resulting  integrals,  and  using 
equation  (28),^  yield  the  Hertz  potentials,  29-38,  which  are 
the  solutions  to  equations  (1-6)  and  (11-14). 


JlpltaL 

V\2-k2 


T^k2 


\dV= 


ik  V r2+z2 


4^? 
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Vertical  Electric  Dipoles 
Case  A 

29.  TT^rjZ,*')  = p 

(A) 


i<».)  t - 

4tts  |_  R * 

rcoshmdj 


2. 


ipa 


t « 


R*i 


^(z+z»-2d)  xJo(*j.)dx, 

d) 


00 


30.  TT2(r,z,zf)  « p^z*)  j cosh 
(A)  „ 2 Jo  . 


JB2-S. 


Region  1 
--£(z  r-d) 


\Jo(*.r)dK 


2nen 


M(0,\,n,d) 


Region  2 


Case  B 

31. 


32. 


Case  C 

33* 


n1(r,z,zl)  *=  n2(r,zf,z)  p2(z*)  = 

(B)  <A) 


<*  7 

en2  Jo 


*2 
2nen2 


cosh  mz»e 
H(p,X,n,d) 


-l(z-d) 


\Jo(X,r)d>. 
Region  1 


lompR  in3R 

Tr2(r,z,zf)  = p2(z» ) le  + e 

(B) 


2 

4rm  e 


R 


- 2 


/&  , -md 

I (‘■t-  ^lyn2)coshI^z,  coshmze  \Jo(\.r)dX. 

W tr/O  a ^ ^ \ 


m M(3,\,n,d) 


Region  2 


n1(r,z,0)  - p2(0) 


(C) 


4rren 


1 


r^(z"d^.jo(\r)d7l 


M(3,X,n,d) 


Region  1 
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34 


. it2(p,z,0)  = p2(0)  einPR  - /<*-“/ n2)coshmz  e“mdKJo(Kr,dX 

4fm2e  R o m M(0,X,n,d) 


Region  2 


Vertical  Magnetic  Dipoles 
Case  A 

35.  ^(r^z*)  = nm^z') 

(A)  4rr_ 


eiPR»  eiPH» 

R‘  R" 


+ 2 


36.  TTni2(r’Z’z:,)  = m 

(A) 


Case  B 

37.  * n^r^'zjm^z') 

(B)  (A) 


7*  sinh  md  e ‘"^Z+Z,”2<3^XJ  (X.r)dX.l 
/ TO^,d)  ° 

Region  1 

,(z')  f slnhmz  e-^z'-d^j^JdA 
2tt  *4  N(P,X.,nd) 


Region  2 


= / slnhmz * e ”^z“^^AJ0(Xr)d\ 

2tt  *4  N<p,\,n,d) 


38.  TT^r^z*)  = *im2(z 


4tt 


2 p1BRI 


Region  1 


- e£fi 

RV 


- 2 


I <fc* 


)sinhmzt slnhmz  e ”mdXJo ( Ar ) dX 


mNO,A,n,d) 


] 


Region  2 


Where 


39.  H(9,X,n,d)  = tcosh  md  + sinh  rad 

40.  N(9,A,n,d)  =^sinh  md  + m ^osh  md. 


s 
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R Is  the  spherical  coordinate  distance  of  Figs.  2a  and  2c. 

=yr777 

R«  ="/  r2  + (z-z«)2  Fig 

R"  =Y  r2  + (z+z’-2d)2  Fig 

RIII * V  =l/  r2  + (z+z»)2 

It  is  to  be  observed  in  equations  29-3^  that  the  forms  of 
the  integrands  of  the  unintegrated  ' terms  indicate  that  these 
are  standing  wave  distributions  or  multiple  reflections  of  the 
waves  from  the  dielectric  and  metal  surface  in  region  2,  the 
dielectric  region.  The  integrated  terms  indicate  that  the 
dipoles  are  Imaged  In  the  dielectric  or  metal  surface. 

The  E and  B fields  may  be  derived  through  the  application 
of  equations  8 and  9 or  16  and  17« 


2a 

2a 


III 


Far-Zone  Formulations  of  the  Hertz  Potentials 

The  Hertz  potentials  and  for  vertical  electric 
dipoles  are  asymptotically  integrated  to  obtain  far-zone 
expressions  with  the  restriction  of  thin  dielectric  coatings. 
The  Integrations  of  equations  2.29-2.31  and  equation  2*33, 

by 


Tai10  and 


Appendix  A,  are  similar  to  those  performed 
11  12 

others.  ’ Equations  2.32  and  2.34  are  integrated  using 
the  method  of  Van  Der  Waerden,^  Appendix  B. 

The  integrals  in  equations  2.29-2.34  are  transformed 
to  a symmetric  form  for  r>  o: 


(1) 


I 


F(\)J0(Xr)M\ 


■if 

'■'—CD 


F(X)H0A(Xr)\d\ 
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through  the  use  of: 


JQ(Xr)  = 2 Cff01(X,r)  + H02c^r):i 


ff02(Xj)  = - H^C-Xj) 


where : 


F(X.)  is  an  even  function  of  X.,  and  H^CXj*)  is  chosen  in 
order  to  yield  outward  travelling  waves  for  real  positive  X.. 

The  integrands  in  equations  2.29  - 2.34  have  first-order 
poles  where  M(0,\,n,d)  = 0,  equation  2.39,  or  for  solutions  of: 


vLp2  = ± 


P2P 

n*0  -X  tan(d 


Examining  equation  (2)  in  the  light  of  equations  2.29  - 2.31, 
it  is  found  that  i and  m must  be  either  positive  real  or  negative 
imaginary  numbers  in  order  not  to  violate  the  radiation  condition.'*’ ^ 
It  Is  also  observed  that  solutions  exist  only  for  the  same  sign 
on  both  sides  of  equation  2 and  that  the  roots  must  satisfy: 

(3)  f>  - Nr*1* 

Figure  4 contains  curves  of  XT/0  versus  0d  the  electrical 

2 " 

dielectric  thickness,  for  n = 2.54,  polystyrene,  and  for 
n2  = 2.25,  polystyrene,  in  the  range  .lOnr<0d  S 1.6n.  Equation 
2 has  at  least  one  solution  for  all  0d,  0d  > 0.  A new  solution 
or  root  X^  appears  when: 

IciY 

(4)  w =-yT2^  k = °,1,2»  • • • a • 

Henceforth  k will  designate  the  number  of  solutions  of 
equation  2 for  a given  0d. 

The  integrands  in  equations  2.29  - 2.31  and  2.33  have 
branch  points  of  i = 0 or  X " + 0,  (1)  and  (2)  in  Fig.  3, 
and  at  X = 0,  (5)  in  Fig.  3,  and  poles  at  X = A^,  . . . , X^. 

Upon  examination  of  the  integrals  for  slightly  complex  0 and  n, 


i 

f 
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It  is  seen  that  the  integrals  are  to  be  taken  over  path  in 
the  complex  X.  plane  of  Fig..  3» 


The  integrands  in  equations  2.32  and  2.34-  have  additional 
branch  points  at  m - 0,  \ = + n{3,  (3)  and  (4)  in  Fig®  39  and  in 
addition  the  same  poles  and  branch  points  as  the  former  integrands. 
These  integrals  are  to  be  taken  over  W2« 

The  evaluation  of  the  k residues  and  the  combination  of  these 
with  the  assymptotic  integrations,  performed  in  Appendices  A and  B, 
yield  the  far-zone  Hertz  potentials,  equations  5 - 10,  subject  tog 


where  R = y r2  + z2 

<<  0(R“2)  <<  OCR"1) 
e“R  |n£  - **|«0(R"2)  ^<0(R“1) 

or  in  other  words  the  poles  are  not  too  near  the  branch  points. 

The  terms  in  the  summations  in  equations  5-10  which  are 
due  to  the  evaluation  of  the  residues,  are  waves  which  propagate 
radially  and  attenuate  exponentially  in  the  z direction.  These 
waves  "hug"  the  surface  and  are  hence  called  guided  or  surface 
waves.  They  result  as  eigenfunctions  of  the  configuration,  Fig.  I, 
(excluding  r = 0).4"*  ^ The  remaining  terms  will  be  designated 
as  the  radiating  or  compensating  waves. 

The  far-zone  Hertz  potentials  ares 


and 


0R» 


1 

and 

Pd 


rRlp  -\\ 


Case  A Hertz  dipole  above  the  dielectric  at  diz'Soq  r'  = 0 


v PiU') 

5«  ni  (r,z,z')  = - 


(A) 


e 


10R 


e 


i0R 


+ KSj).  e 


4ttR  4rrR 


1/Z^jT 


d < z H oo 


1 i v 


FIG.  3 PATHS  OF  INTEGRATION  IN  THE  X=  ij  + ts  PLANE 
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P P-.Cz') 

6.  n2F(r,z,z»)  = — y — 1 

(A)  n e 


P(©2)cos(pzVn2-cos2©2)  W’ 
= — e 

R,,f  cos(pd  Vn2-cos2©2) 


Aj(e2)Q(\J$3,n9d)cos(zV  n2p2-*.j2)  CiV~(z'“d)  V 


j=1  V2ttA.jP  cos  (d  V n p “A-j^ ) 

0 < z <d 

Case  B Hertz  dipole  in  the  dielectric  at  0 ^ z<d,  r'  = 0 


„ P0(z!)  P(©,)  cos  ( pz » V n2-cos2©-. ) i0R 

7_®  nnF (r,z,z* ) = — p — 1 3 / =~  e 

(Pd  V n2-cos20-j) 


R^  cos 


+ ^y7  AJ(Q^)Q(XJ,p,n,d)cos(z»  V n2p2-*-j2)  CiA-jr-Cz-d)  ]/ \j2-p2] 
V2tt\jP  cos(d  / n^P^-Aj^) 


d £ z < oo 


w PP(z‘) 

8,  it2  (r,z,z‘)  * ■ — 

(B) 


en 


QQ^P^dJcosCz1  *[/"n^p<::-A.j£:)co3(z y n^p^-xy) 


.2a2  , 2, 


l/n^p2-*..,.2)  iX,Jr 


=£  l^nXjT 


cos 


5(d V n2p2-\2) 


e 


'J 

0 2S  z <d 


Case  C Abratiam  dipole  on  the  metal  surface  at  z’  = 0,  r* 

,IV 


= 0. 


P P2{0)  , 

9.  ny<r,z,0)  = 


(C) 


iPR 


Aj(©^)Q(Xj9P9n9d) 


P(Q^)  e 

R1^  cos(Pd  Vn^-cos2©^) 

(iXTr-(z-d)  lOp2)] 


2n\^  cos(d"/  n2p2-\j2) 


d < z is  oo 
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» P o<0)  Qaj,3,n,d)cos(zl/  n2?2~\  2)  iLr 

10 o TTg  (r^z^O)  - ^ ® 

(6  2n  e j=l  Y 2n^r  oos2(a  y ^V^2) 

0 < z <d 


P(e)  = 


i sin€ 


2n[i  sin©  + n"2  V n2-cos2©  tan(  3d  V n2-cos2©)  ] 

-iTT/4  "^.j2-p2.  (n2p2-^j2) 


Q(^j9  P^n^d ) 


02(n2-l)  + d lA.T232[n2(\T2-32)+( p2~\  2n"2) ] 


and  Aj  (0)  is  defined  byj 


Aj(6) 

= 0 

© > 

cos 

_1PAj 

Ve) 

...  1 
2 

© = 

cos 

_1PA, 

Aj(8) 

= 1 

© < 

cos 

_1PAj 

sin"1 

iiLzSjU 

R 

s cos-1  r/R, , 

»»  t 

R = T 

^r2+(z+z 1 

-2d)2 

sin"1^,^ 

R 

= cos"1 

r/R  ‘ ' 

< « 
R 

i _ • 

j/  r2+(z' 

-a)2 

0,  = sin"1  = cos"1  - 

3 R*lv  R 

Rw  = Vr2+(z-d)2 


Pig.  2a 


Pig.  2b 
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« * * 

If  it  is  assumed  that  0z'<<0R  , the  first  term  in  6 

vanishes  and  Aj(©2)  = 1. 

The  magnitudes  of  the  fields,  of  the  guided  modes  attenuate 
with  height  z in  two  ways  - first,  the  exponential  or  guiding 
factor  exp(-z"|/  \j2-02)$  second,  the  coefficient  Aj(6).  The 
second  factor  is  a function  of  0r  while  the  first  factor  does  not 
depend  upon  the  radius.  It  may  be  observed  that  as  0r  increases  the 
exponential  attenuation  takes  precedence  over  the  coefficient 1 s 
attenuation.  Pig.  5 contains  curves  of  the  minimum  radii  0r, 
for  each  mode,  at  which  the  field  attenuates  to  1 percent  of 
Its  maximum  value  due  to  exponential  attenuation,  at  a height 
0z  less  than  that  at”  which  the  coefficient  Aj(0)  equals  its 
critical  value.  These  0r  are  plotted  versus  0d,  in  the  range 

p T 

. IOtt  ^ 0d  ^ 1.6tt  and  for  n = 2.54  and  2.25.  For  0r  greater 
than  these  values  Aj(©)  s 1 to  within  1 percent  accuracy. 

The  presence  of  the  factor  Aj(0)  indicates  that  as  0r 
increases  the  guided  wave  acquires  more  power  until  at  0r  * co 
the  total  power  alloted  to  it  has  been  attained.  Making  the 
assumption  that  Aj(0)  Is  identically  one,  for  0r  greater  than 
the  values  of  Pig.  5,  Is  the  same  as  assuming  that  the  guided 
wave  has  acquired  its  total  power  at  this  0r.  It  is  valid  to 
within  .01  percent  accuracy  with  respect  to  the  Poynting  vector. 

Since  in  this  report  the  far  fields  are  being  investigated,' 
it  will  be  assumed  henceforth  that  0r(=0R  sin©)  is  greater 
than  the  values  of  Pig.  5 for  each  0d  being  investigated.  For 
this  reason  the  notation  Aj(9)  will  be  dropped  under  the 
assumption  that  Aj(9)  = 1. 

It  is  to  be  observed  that,  in  order  for  the  fields  to 

p 

be  identical  for  sources  at  z*  = d+  and  z*  = d-,’  p^(d+)=p2(d-)/n  . 
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IV 

The  Far -Zone  E-  and  B-  Fields 

The  representation  of  the  far-zone  E-  and  B-  fields  are 
obtained  by  applying  equations  2.8  and  2.9  to  equations  3.5  - 
3ol0o  Terms  of  order  higher  than  r“^  are  neglected,,  The 
assumption  is  made  that  0z'<<  pR  so  that  the  radiation  term 
in  equation  3°6  may  be  neglected,, 

It  is  convenient  to  represent  the  radiating  terms  in 
spherical  coordinates „ The  following  substitutions^  are  made 
in  the  radiating  waves  assuming  0R»1  in  Fig.  2. 

In  the  phase  factors: 

0R*  = 0R  r Pz'  cos®  pR"  = 0R  - p(2d-z*)cos0 

TV1 

PR-*-*  * pR  - pdcos© 

In  the  amplitude  factors: 

r/R"  = cosOj  = sin©  = r/R  = r/R* 

ft"  ■f|  n ' = slnO^  « COS©  = z/R  = (z  - 2*^1  * 

PR'  = PR"  = pR 

sinO^  = r/R17  ^ r/R  = cos© 

cos©2  = r/R*V  = r/R  = sin© 

PR  = pR17 

The  substitution  of  these  into  equations  3*5  - 3.10  after 
2.8  and  2.9  have  been  applied  leads  to  the  following  far-zone 
fields  for  vertical  dipoles.  These  are  subject  to  the  same 
restrictions  as  3.5  - 3.10  in  addition  to  Pz'<<ipB. 

Case  A Hertz  dipole  above  the  dielectric  at 
d^  z*  — od,  r*  = 0. 


FIG.  5 MINIMUM  RADIAL  DISTANCE  $r  AT  WHICH  THE  Au  (0) 
COEFFICIENTS  OF  THE  SURFACE  MODES  ARE  IDENTICALLY 
1 FOR  ALL  Z AS  A FUNCTION  OF  DIELECTRIC  THICKNESS 
$6  FOR  0<  /3d  < 1.6  TT.  BASED  UPON  AN  ACCURACY  OF 
1%  OF  THE  FIELDS  AND  0.01%  OF  THE  POYNTING'S  VECTOR 
n2  - 2.54 , n2  = 2.25 
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Region  1. 


1. 


.A 

.R  x B1  . 
(A) 


2. 


_*F 

_J?r 

k 

E^r^z1 

') 

= Ei 

+ )>  v 

E1J 

(A) 

(A) 

(A) 

Ci 

— 5 -»FG 

-r  Bu 

x £ 

— FG 

- r x Aj  B:j] 

0 J=1 

(A) 

(A) 

F 

— Fr 

Bl(r,z,z 

') 

= B1 

+ A 

Bu  ‘ 

CA) 

(A) 

J=1 

(A) 

a p^Cz' )3^sln0  - CcosBcos^  + S sinjf]  10(R-dcos8] 
i 2ttRu)e  f cos©  - IS  1 6 


,2  k 


© p,  (z'  )0  CiXT 

-iy^5w,M,a)e 


Cik.Tr-(z+z,-2d)  i-K  t2-32] 


Where:  R,  ©,  and  $ are  spherical  coordinates „ 

r,  ©,  and  £ are  cylindrical  coordinates. 

_ 1 

V * (eii)“  2 = the  characteristic  >rslocity  of  medium  I 
- ^(B,A,z,,6)  = PCd-z'Jcos©. 

S * S(n,0,d,0)  = n“2*l/  n2-sin2e  tan(0d"|/  n2-sin20). 


Region  2, 


3. 


_*F 

B2(r,z,z') 
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E2J 
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= ~~22 
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i V~ti 


202-Aj2  tanC  z V n202-\T2)B 
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P \ / t?G 

B'2(r,z,z»)  = B2J(r,z,z«) 

(A)  J=i 


©P-.Cz1)^  \~~7  ‘■-‘-"•T 

= - — 1 . > \TQ(X  3,n,d)e  J 

we  fzf 


dX.Tr-(z!-d)  y X.t2-32]  • 


Ccos(z  ~\f  n2g2-X.j2] 


CcosCd"/  n2p2-k  2] 


Case  B-  Hertz  dipole  in  the  dielectric  at  0<  z*£d,  r' 
Region  1. 

_^F  _^Fr  v ^ / 

5»  E1(r,z,z ’ ) -E1  (r,z,z«)  + p>  ^ Eu(.r,z,z*) 

(B)  (B)  J5!  (B) 

«.VT  r ^ 7 / 7*j  7Z  _JG  . . — JG 

= - V(R  x B1  ) + JT  ClV  >.j2-0  Bu-X  z - r x V,  B1;r) 
CB)  P W (B)  (B) 


B1(r,z,z‘ ) - Bx  (r,z,z')  + B^j(r,z,z«) 

(B)  (E)  (B) 


_ _ ^ p2^z'  ^ CslnBcosQ][cos(Bzl  “]f  n2-sin26)  dcosQ)j 

2Twn2eR  [cog0.iS]  CcosOdV  n^sin2©)  ] 


A 

- e 


S • 


[ iXjr- ( z -d  )Yx./-P2]  • 


Region  2 

-pr 

E2(r,z,z») 

(B) 


cos(z‘ V n2B2-Xj2 ) 

cos(d  V n202-XT5) 


^F6 

B2J(r»z’z,) 
W (B) 
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= V;  CiV  n202-\j^  tan(zV  n202-X.j2)  B2J  x z - $ x X.jB2j3 

Tl  ft  • ■ ^ # v ^ _ \ 


8, 


* > -J=T 
F 


(B) 


F FG 

B2(r,z,z«)  = > B2j( 

(B)  “ (B) 


r,z,z‘ ) 


*P2(z')02 
am  e / 2ttt  J=1 


d)e 


i\jT 


cos( z * V n2p2"\j2)cos(z  V n2p2~X.j2) 


cos 


2(dl4232-X.T2) 


(B) 


Case  C.  Abraham  dipole  on  the  metal  at  z ’ = 0 , r1  =0, 

The  respective  quantities  are  given  by  the  following 
relations; 


_JF  i 

9-  B (r,z,0)  = f 

(C)  2 


lim  E (r.z.z1) 
- <B)  ’ 

Z4  0 


10 


. ^ (r,z,0)  = 4 lim  ^ (r  z 

(C)  2 z'-^o  (B) 

with  subscripts  1 or  2. 


z* ) 


It  is  to  be  noted  that  the  radiating  or  compensating  terms 
are  waves  of  spherical  type  which  attenuate  as  the  guided 

terms  are  waves  of  cylindrical  type  which  attenuate  as  r"^.  As 
r or  R increases,  the  ratio  of  the  magnitude  of  the  guidea  wave 
to  that  of  the  radiated  wave  increases. 


The  guided  modes  are  elliptically  polarized  in  region  1 with 
semi-majjor  axis  in  the  z direction  and  semi-minor  axis  along  r. 
The  eccentricity  for  each  mode,  which  is  Independent  of  position, 
is: 


11. 


e - 3Aj 
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The  guided  modes  are  sometimes  referred  to  as  "slow  waves" 
because  their  phase  velocities,  Vj  = w/X.j,  are  less  than  the 
characteristic  velocities  of  the  region  1. 

Goubau^’1'7  and  Attwood4,  have  treated,  at  length,  the 
field  and  flux  distributions  for  waves  of  this  type. 


Field  Patterns 

The  field  patterns  are  curves  of  the  relative  magnitudes 
(with  respect  to  their  maxima)  of  the  6 -components  of  the  E- 
fields  in  region  1,  as  functions  of  g,  in  the  range  of  0 between 
0 and  tt/2,  for  0R>>  1.  The  dielectric  is  assumed  to  be  thin 
so  that  0 = tt/2  along  the  surface  between  region  1 and  2,  Fig.  2* 
To  obtain  the  formulas  from  which  these  curves  are  computed,  it 
is  necessary  to  convert  the  guided  components  of  the  fields, 
equations  4.1  - 4.10  to  spherical  coordinates  and  to  recombine 
the  fields  at  some  large  fixed  0R. 

It  is  important  to  notice  that  the  field  patterns  defined 
above  are  those  measured  by  a receiving  antenna  polarized  and 
traveling  in  the  6-direction  at  the  chosen  0R.  Since  the  guided 
waves  have  R components,  the  patterns  are  not  the  relative 
magnitudes  of  the  total  field  along  the  great  circle  of  radius 
0R.  It  isalso  important  to  observe  the  manner  in  which  these 
patterns  vary  with  0R.  Since  the  radiated  wave  disappears  at 
0 = tt/2,  while  for  most  finite  dielectric  thicknesses,  0d,  the 
surface  wave  attenuates  rapidly  with  decreasing  0,  the  region 
of  interference  between  the  two  is  very  small.  For  these 
reasons,  the  principal  effect  of  varying  0R,  is  to  vary  the 
relative  magnitudes  of  the  sections  of  the  field  patterns,  at 
0 near  tt/2  and  at  0 less  than  tt/2,'  by  a factor  R^. 

Figures  6-14  contain  field  patterns  for  cases  A,  B, 
and  C,  n2  = 2.54,  0R  = 239  and  various  0d  and  0z’,  on  semi- 
logarithmic  scales,  with  comparison  patterns  for  dipoles  above 
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FIG.  6 FIELD  PATTERNS  FOR  CASE  A,  = O.IOTT,  M(O)  VS.0 


WITH  FIELD  PATTERNS  WITH  £d»0  FOR  COMPARISON. 

£R0  = 239,  n2*  2.54 , £d  = O.IO  T ,/3d  = 0 

/3z'  * O.IOTT,  TT/4  +•  O.IOTT  ,TT/2  + 0. 10 TT,T  + O.IOTT 


FIG.  7 FIELD  PATTERNS  FOR  CASE  A,  £d  = O.l 5*JT,  M (0)  VS.  0 
WITH  FIELD  PATTERNS  WITH  £d  = 0 FOR  COMPARISON 

/9R0=  239,  n2  = 2.54,  £d  = 0.I5T , 0d  = 0 

/3z'»  o.i5*rr,TT/4  +o.i5T,ir/2+o.i5'ir,,ir+o.i5Tr 


FIG.  8 FIELD  PATTERNS  FOR  CASE  A,  £d  = nT/4,M(0)  VS.  0, 
WITH  FIELD  PATTERNS  WITH  fid  = 0 FOR  COMPARISON 

£R0=239,  n2=2.54,£d  = nT/4 , fid  - 0 

fiz'  = nT/4  ,(n+ 1 )T/4 , ( n + 2 ) T/4 , ( n + 4 ) T/4 


FIG.  9 FIELD  PATTERNS  FOR  CASE  A , /3d  WITH 

FIELD  PATTERNS  WITH  /3d  =0  FOR  COMPARISON 
/3R0  = 239  ,n2=254,  /3d  = ^-,  /8d  = 0-—  , 

/3z'=^,-^+^,(n+l)-|,  (n  + 2)-f- 


I 


FIG-  10  FIELD  PATTERNS  FOR  CASE  B,  /3d  = O.IOTT,  n(0)  VS.  0 
WITH  HELD  PATTERNS'  WITH  /3d=0  FOR  COMPARISON. 

/9R0=  239, n2=  2.54,/3d  = O.IOTr ,/3d  = 0 

$z'  * 0.  IOTT/4 , 0.  IOT/2 , 0. 10  x 3T/4 


FIG.  II  FIELD  PATTERNS  FOR  CASE  B,  £d=O.I5ir,  M (0)  VS  0 

WITH  FIELD  PATTERNS  WITH  /9d  = 0 FOR  COMPARISON 

£R  =239,  n = 2.54, 0d  = 0.15*—,  £d  = 0 

£z'  = 0.I5-J-,  O.I5-|-,O.I5x^ 


j8R0=  239,  n2=254,  = £d  = 0 — 


$H  • 239, n =2.54,  /3d  =0 — ,£d=O.IOir,O.I5T, — 
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or  on  a perfect,  uncoated  conductor  for  the  same  dipole  heights, 
Pz’,  as  those  with  the  coating  present e 

The  dielectric  constant  of  polystyrene  was  chosen  as  being 
typical.  0R  = 239  was  not  only  chosen  as  being  typical  but  was 
chosen  for  proposed  future  experimental  use.  The  thickness 
Pd  = nrr/2  may  seem  peculiar,  but  in  Pig.  4 this  value  is  seen 
to  be  smaller  than,  but  very  close  to,  tt/V^  the  value  at 

which  multimode  propagation  commences.  Since  the  value  Pd  = nrr/2 
possesses  this  quality,  besides  being  convenient  for  numerical 
computation,  it  was  chosen.  It  is  clear  that  the  dielectric 
thicknesses  used  provide  good  coverage  of  the  region  of  single- 
mode propagation. 

The  quantities  plotted  in  Pigs.  6-14  ares 

I Fl 

1.  M(6 ) = Ee  | 

|®oJ  0 ~ n/2  ~ ®max 
For  Pd  = 0 and  Pd  f 0. 

The  main  points  of  interest  in  Pigs.  6-14  are  the  large 
maxima  or  spikes  at  0 5 rr/2  and  the  compressions  of  the  rest  of 
the  fields  with  the  accompanying  minima  between  these  two.  These 
curves  are  drawn  on  a logarithmic  scale  and  must  be  examined 
closely  to  appreciate  the  actual  large  magnitudes  of  the  maxima 
at  0 = tr/2  (the  surface  waves)  compared  to  the  rest  of  the  pat- 
terns (the  radiated  waves).  The  ratios  of  the  powers  in  the  two 
components  of  the  fields,  Pigs.  20  and  21,  give  a better  concep- 
tion of  the  large  size  of  these  maxima,  especially  if  consider- 
ation is  given  to  the  fact  that  the  surface  wave  is  compressed 
into  a small  segment  of  space,  while  the  radiated  wave  is  dis- 
persed over  the  complete  half -space. 

For  case  A the  field  patterns  for  0 < tt/2  are  similar  to 
those  without  the  dielectric  coating  for  dipole  heights  equal  to 
the  dipole  heights  in  case  A minus  the  dielectric  thickness. 
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This  means  that,  for  case  A,  the  radiated  waves  are  strongly 
reflected  by  the  dielectric  surface. 


In  cases  B and  C the  effect  of  the  dielectric  is  to  compress 
the  radiated  wave  into  a smaller  sector. 


Figure  15  is  a plot  of |Eq 
Figs.  6 


FQ 


14  over  the  range  of  significant  magnitudes  of  ■5° 


versus  G for  the 


Pd's  of 
> For 


the  cases  considered  S’  - 1.  These  yield  an  indication  of  the  magni- 
tudes of  the  R-components  of  the  radiation  patterns.  For  the  thicker 
dielectrics  these  ratios  are  almost  constant,  indicating  that,  since 

Fig, 

tively  over  the  range  of  significant  values  of 


/ E~, t is  constant  and  equal  to  the  0 to  R component  ratios  of 

/ « T ?P  W1  wp 

15  at  G * n/2,’  E^j  and  E^J.  are  essentially  and  E*  j respec- 


VI 


Formulation  of  the  Power  Equation 

The  purpose  of  this  section  is  to  formulate  the  equations 
which  govern  the  time-average  flow  of  power  across  a surface 
surrounding  the  source  p(z').  Since  the  far-zone  approximations 
do  not  satisfy  Maxwell's  equations  and  since  the  proofs  that  the 
powers  in  the  guided-and  radiated-type  waves  are  independent  re- 
quire the  fields  to  satisfy  Maxwell's  equations,  it  is  necessary 
to  formulate  the  exact  fields. 

Consider  E andHT,'  the  exact  fields  existing  in  the  configu- 
ration of  Fig.  1,'  due  to  a source  p(z*)  in  any  of  three  positions 
A,  B,  or  C of  Fig.  1. 

For  0r>  0 the  transformation  of  equations  2.29  - 2.34  to  the 
Hankel  function  form  3.1,  is  valid.  The  exact  integrations.  If 
they  could  be  performed,  would  be  taken  over  paths  and  Wg  of 
Fig.  3.  These  integrations  yield  contributions  from  the  poles  and 
branch  cuts.  The  application  of  equations  2.8  and  2.9  to  the  re- 
sults of  the  integrations  yields  the  exact  F-and  B-fields.  It  is 
sufficient  for  the  purposes  of  this  section  to  represent  the  branch- 
cut  contributions.  Er,~Sr,  symbolically,  while  the  contributions 
of  the  poles, "3^,  may  be  explicitly  evaluated. 


FIG.  15  Eq/Er  VS.  0 OVER  THE  RANGE  OF  SIGNIFICANT 

values  of|eg|  for  £d  * nT/2,  mr/4,o.i5Tr,o.io*ir 

AND  FOR  n2=  2.54 


1 
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Por  pr>  0 the  exact  fields  may  be  represented  bys 


1.  a)  *E  =~Er  +~iG 

2. 


b)  B = F + Bl 


^ C 1? 


J = 1 


j=i 


/\  / o ? 1 a.  t -z^C^-p 

rr-L/-\  -\  x rr-L/\  J 


2 «2 


Au  ® [r  hJ(v)  + zKjHo(xJr)]e 


b)  Bl 


B° 

BJ 


J=1 


i «i 


B2  J 

<v>  t 


J=1 


Region  1 

i2S2-\j2  sin(zl/n2B2-X.j2 

HiX(V> 

Ho<V>] 

Region  2 

-z  Jx.  2-p2 

A ■ V 

U ® 

Region  1 

2j  n2cos(z  Vn232-^»j2) 

Region  2 

Where  Aj  is  a function  of  X.J9  B,  n,  d,  a>9  z»  and  p(z‘). 
For  3B>>1 
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Where  the  fields  on  the  right  sides  are  given  in  section  4. 

and  are  approximate  solutions  to  Maxwell's  equations 
for  large  0r  while  E,  B,  Er,  B1,  E^,  Bj  are  exact  solutions 
to  Maxwell's  equations  for  3R>0. 

For  0r»l  the  guided  or  residue  waves  exist  uncancelled, 
hut  as  0r-*-O  the  functional  forms  of  the  guided  waves  become 
infinite.  Since  this  is  physically  impossible  except  at  discrete 
points  at  r = o,  z = z',  the  so-eAlled  radiating  or  compensating 
fields,  due  to  the  branch  cut  contributions,  must  cancel  the 
singular  parts  of  the  guided  waves  as  $r— ►O.  An  indication  cf 
the  cancellation  is  given  by  the  appearance  of  the  factor  Aj(0) 
in  the  asymptotic  results.  Despite  this  cancellation  the  division 
in  equation  1 is  valid  because  of  the  linearity  of  Maxwell's 
equations. 

Since  p(z')  is  arbitrary,  let  it  be  chosen  to  be  a real 
quantity.  It  may  be  verified  that  for  p(z')  real  A-^j  and  A2j 
are  real. 

The  power  flowing  across  any  surface  S surrounding 
the  source  is  given  bys 

n c vEr  x 1§p*)dS 


where  ft  is  the  outward  normal  to  the  volume  enclosed  by  S, 
and  superscript*  indicates  the  complex  conjugate. 

It  will  be  proved  that 

5.  §|  C /ft  . (Er  x x~Br*)ds]  = 0 
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(ITj0  x'Bm°*)dS  = 0 


(The  author  is  indebted  to  Dr.  George  Goubau  for 
suggesting  these  proofs  and  recommending  references 
to  similar  proofs.1®’  1^). 
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Some  necessary  formulas  derived  by  well-known  methods  are 


% / * B?dS  ' */V  * 


/* 


S V 

= dV 

where  S is  a closed  surface  surrounding  a volume  V. 

EA9  Ba  and  Eg,  Bg  are  two  independent  solutions 

to  Maxwell's  equation  in  V. 

P and  J are  the  impressed  polarization  and 

current  densities. 

e and  \i  are  assumed  to  be  real  in  V. 

n is  the  outward  normal  to  V. 

21 

and  the  reciprocity  condition. 

• [®A  X*B  X_®A]  d S = 


- top.  / ( . Eg  - Pg  . EA)dV 


JB  * EA  " JA  • ®B^dV 


Equation  6 will  be  proved  firsts 

To  this  end  consider  the  configuration  of  Pig.  16  agreeing 
with  Pig.  1. . 

Since  the  cylinders  S^  and  S2  are  semi-infinite  and  the 
surface  z - 0 is  a perfect  conductor,  S may  be  set  equal  to  S^ 
or  S2  in  equation  (4).  It  is  therefore  sufficient  to  prove 


-30- 
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equation  (6)  over  S1  or  S2* 

If  = Ej^  and  = ~B^  are  substituted  in  equation  (7) 
with  the  volume  as  the  enclosed  volume  and  noting  that  the 
volume  V^,  Fig*  16,  is  3ourceless,  the  following  identity  is 
arrived  at* 


Se  f 
2\i  J 


a ,— *■  0 — »G  .* 

r * (Ej  x Bffi  ;dS1 


% f 


r • (Ej®  x Bm^*)dS2  = constant  in  r = K. 


This  is  a constant  since  r^  and  r2  may  be  independently 
varied  while  the  other  is  fixed  or  explicitly 


10, 


Re 

25 


i- 


r . dT,G  xlTj^dS 


’m  '““1  or  2 


S1  or  S2 


- Be  ” [1  ra  or  2 Vo^Vl  or  2)Hl2(7mrl  or  2)]DJm 


= K 


where  the  constant  Djin  is 
oo 

uJm  to 


J AU  Alm  expC-z(  /Xj2-02  + /xm2-p2)]dz 

ra  , i 

+ j n2A2jA2mcos(z  n232-Xj2)cos(zVn202-\m2)dz 

A ' 


Specifically,  this  is  true  for  large  0r  where  the  asymptotic 
forms  of  the  Hankel  functions  are  valid: 


11* 


lim  K - DTmJ cos(r(\T-V_)) 
pr»l 


For  \j  f this  is  zero  for  r 


qnr 


where  q = integers,  but  K is  a constant;  hence  the  constant 
must  be  zero.  Thus  equation  (6)  is  verified,  since  Djm,  J f m, 
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must  be  zero  for  K to  be  zero  for  all  r > o 


Equation  (5)  will  now  be  proved: 

Consider  the  configuration  of  Fig.  17«  The  source  p(z') 
may  be  in  region  1 or  2 or  half-submerged  in  the  conductor  at 
z = o,  i.e.,  an  Abraham  dipole.  The  current  I(z")  on  the  semi- 
infinite cylinder  V2  is  such  as  to  excite  only  the  Jth  guided  mode, 

After  utilizing  the  representations  of  the  Hankel  functions 
for  small  arguments  the  following  representation  for  I(z")  is 
arrived  at.  / — 5 — 5 

a e -z"  ™-r  -e 

2 pje 


12-  - - iQr 


p 

A2j  n cos 


(z“  Vr?32-A.T2) 


d zn  < 00  . 
6 iz"  < d. 


— ^ p -r  n 

This  current  excites  only  Ej  (r1?z)  and  Bj  (r^).  The  dipole 
excites  E(rQ,z)  and  B(rQ,z). 

The  subscripts  for  cases  A,  B,  or  C agree  throughout.  It  is 
assumed  that  Iz(z")  is  the  total  current  on  F2  and  p(z')  is  the 
total  polarization  in  V^. 

The  application  of  equation  (8)  to  V^,  F^,  noting  that  the  - 
polarization,  is  a constant  times  a delta-function,  leads  to 


13. 


/*' 

F1 


[E  x BjG  - 1®  x B]dF1  = 


cop.  X.jHq1(X.jX) 


-z’//\  2-{32 

fP^z*  JA^j  e 

p p(z* )A2jCos(z 1 Vn232-^j^) 
P2(0)A2j 


d £z'<-oo 
0 ^.z’<.d 

z'=  0 


The  application  of  equation  (8)  to  V2,  F2  in  the  light  of 
the  orthogonality  condition,  equation  6,  gives 
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o CE  x Bj  - x B]d  F2  = 


H J I(z« 


)E  r(x,z“ )dz” 


^2.  2 . . _2  A2J2sln(2d  Vn2g2  - Nr,2) 


- Ho1<V)  n2A2J2  d + * 


2Vn232-\j2 


A,  t2  e-  2d  At2  " P2 


\j2  “ 0‘ 


The  application  of  equation  (8)  to  the  sourceless  region  V^, 
noting  that  the  integrations  over  the  surfaces  at  infinity  and 
at  z - 0 vanish,  and  rearranging  the  resulting  equation  leads  tog 


H / I(z")Ez  dz” 


Ho  <V° 


^ [n2A2d  + n2 


hh  + 


-2dlk  2-02 


P 1<  z ' ) Ajj.  e 


- z'Vy  - r 


d^  z' oo 


+ cop.A.j  I P 2(z')  Agj  cos(z'  Vn202  - Xj^) 


0 < z'<d 


Po(0)A.. 


z'  = 0 


SUPPORTING  CURRENT  Iz(Z")  AND  SPHERICAL 
SURFACE  SURROUNDING  p(Z')  SEPARATED 
BY  DISTANCE  x 
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The  right  side  of  equation  (15)  is  not  a function  x or  z- 
Hence,  the  equation  is  equal  to  a constant,  G,  in  x and  z- 


As  0x  becomes  large,  0x  >>  1$  the  quantities  on  the  left 
side  of  equation  (15)  approach  zero  in  the  following  manners 


1 

/ O 0 

ixd+  K"d 

1 

Vx 

Hence: 

16- 


lim  G = 
0x — >-oo 


lim 

0X ^00 


GO 

[ 


l(z") Ezrdz" 


H 


o<V) 


Since  G is  a constant,  then  G = 0 and  both  sides  of  (15)  are 
identically  zero. 

It  is  enlightening  to  observe  that  the  right  side  of  (15) 
when  .set  equal  to  zero  offers  an  alternate  method  of  determining 
the  Aj-'s.  The  author  has  verified  the  fact  that  the  Aj*s  are 
identical  when  derived  by  both  methods-  This  offers  a check  on 
the  paths  of  integration  chosen  in  Pig-  3® 


Since  the  denominator  of  the  right  side  of  equation  (16)  is 
well  behaved,' AjX>0,  the  numerator  is  zero* 

With  the  substitution  of  the  right  side  of  (12)  for  I(zw) 
the  numerator  on  the  right  side  of  equation  16  becomes g 


17- 


1 


e£(x,zm)cos(z" 


/n202-X.j2)dz" 


0. 


The  following  integral  is  to  be  considered  with  reference 
to  Pig.  16. 
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ft 


°1 

or 

2 
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[E^r^z')  x BjG*(r,z,z')]cl  S1 

or 

2 


Rerr 


r 

* 

£ 
ri  « 

Hi1(\rri) 

‘ ^ 

or 

or 

2 

2 J 

L1J 


oo  -z 


-Z  *j\ 2 ~ 32 


dz 


or  r. 


+ n2A 


2J  j "z 
o 


cos(z  Vn202  -X.2)az 


r = r1  or  r2 


In  the  light  of  equation  (17),  equation  (18)  vanishes. 

Consider  the  sourceless  region  V,  of  Pig.  16  with  the  fields 
E and  B existing  in  the  region.  The  application  of  equation  (7) 
to  the  two  components  Er,  ~Br  and  E^G,'  B^®  of  E and  B with  the 
result  of  equation  (18)  yields. 


19. 


f T . (S'/  X = §*■ 


s. 


f r . (T/  + "B*r)dS2 


S, 


= constant  In  r 


Since  r^  and  r2  may  be  independently  varied. 

p 

The  substitution  of  the  explicit  expression  for  Ej,  allowing  0r 
to  become  very  large  and  examining  the  result  in  a manner  similar 
to  that  used  in  deducing  that  the  constant  inequation  (15)  was 
zero,  verifies  that  equation  (19)  is  zero  for  all  r. 

Combining  the  results  of  equation  equations  (18)  and  (19) 
verifies  equation  (5)« 

It  follows  that  equation  4 becomes 


t 
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p s E®. 

F 2n 


(&  xt^dS 


i { a • * B°*)ds 


- p3-  + p° 

This  allows  the  division  of  the  power  transferred  across 
a surface,  S,  surrounding  the  source  into  two  or  more  parts 
and  hence  to  separately  compute  the  powers  radiated  and  that 
transferred  across  S by  each  guided  mode. 


The  time-average  power  output  of  the  sources  may  be  computed 
conveniently  by  one  of  two  methods: 

a)  the  integration  of  the  Poynting's  vector  over  a surface 

surrounding  the  source  in  the  far  zone:  or  b)  an  integration 
at  the  source*  Both  methods  lead  to  the  same  results.  Method 

b)  is  preferred  for  simplicity. 

Let  use  be  made  of  equation  (6.7)  with  E&  * E,  Bg=B  and 
Pg  * P,  noting  as  in  (6.13)  that  the  polarization,*  P,*  has  the 
form  of  a constant  times  a delta-function  to  obtains 


P = | n . d x~§)dS  = - Re  ^ J 7p*  . 


E av 


= Un  ^ 


p(z')*  Ez(0,z») 
Ez(0,0) 


Hertzian  dipole  at 
r*  = 0,  z*  f 0 

Abraham  dipole 
r ' = 0 * z'  =0 
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Where  S is  a surface  surrounding  the  source  and  V is  the 
volume  enclosed  by  S. 

If  E_(0,z')  is  computed  from  equations  (2*29)  to  (2.34)  by 
the  application  of  equations  (2.8)  and  (2.9),  and  substituted  in 
equation  (1),  and  use  is  made  of  'Hospital's  Rule  for  evaluating 
indeterminate  fractions  or  the  well-known  expression  for  the  power 
radiated  by  an  isolated  Hertz  dipole2^’  24-  is  utilized,  the  ex- 
pressions for  P result: 

Case  A0 


where: 
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If  the  Integrals  are  closely  examined  and  integrated 
wherever  possible  they  are  seen  to  result  in: 

Case  A 

5.  Bn  T I e_2^z,‘"d)dK  = 


x3dx(v^xgcosC2B(z-d)Vl-x2]l-n“Vn2-x2tan(  Bd^-x2)sln[2B(z'-d>ji-x2I 


i“4(n2-x2)tan2( pd\42-x2)  + 1 -x2 


+ n V7- e"ln/4  X.2  e“2(z,-d)  Al*  " ^ Q(X  ; 0,  n,  d) 


7 


J=] 


Case  B 

OD 


6 im  | I »C£m~1  - n~2)e“md  cosh2mz  'dX.  _ 
. r?cosh  md 


l 

^1  _ sln(28z»  , 

Afl*  * t3  Ar3 


4Pz*n^  6n^ 

.+  cos(2flz'^-l)  . sln(2ez'/i?-I) 

4n^02z'2  8ft3*.i3n3 


8pJz,:5nJ 


7 


Vi^? 


1 /-  g3dm 


-r 


-2 


COS2(0zV^) 


n~^i2-x2  tan2  ( 0d/ n2-x2 ) ^/l-x2 


n_4( n^-x2 ) tan2( dVn  -x^ ) + 1-x^ 


$ 


X.j2e"in/4cos2(z'i/n2p2-\  2)Q(\Jf0,n,d) 
■■ - 


cos2(d/n202  - X^) 
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Case  c 


where  Q(A.j,"0,n,d)  is  defined  in  section  III. 


■ „ ■ ■ - < 
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If  method  (a),  in  conjunction  with  equations  (4.20),  (4*1) 
to  (4.10)  were  used  to  compute  PG,  it  would  be  found  that  this 
yields  exactly  the  last  term  in  equations  ( 5“7)  =>  Hence,  the 
terms  in  the  summations  when  substituted  in  equations  (2-4)  are 
the  total  powers  transferred  into  the  guided  modes  from  the  dipoles 
It  is  to  be  noted  that  it  is  necessary  to  prove  that  the  power 
can  be  subdivided  as  in  (4.20)  before  the  two  divisions  can  be 
compared. 

As  was  previously  observed  the  total  power  allotted  to  a 
guided  mode  (the  summations  terms  in  equations  (5-7))  is  not 
acquired  by  the  mode  until  0r  becomes  infinite.  To  a good 
approximation  this  is  true  for  fir  greater  than  the  values  of  Fig.  5 
To  this  approximation,  then,  equations  (2-7)  may  be  given  the 
following  interpretation  according  to  equation  (4.20).  P is  the 
total  power  transferred  across  S. 

J=1  is  equal  to  the  summation 

terms  in  equations  (2-7)  and  is  equal  to  the  power  transferred 
across  S by  the  guided  modes.,  where  S has  a minimum  radius 
0R  which  is  greater  than  the  values  of  0r  in  Fig.  5.  P1*  is  equal 

to  the  remaining  terms  of  equations  (2-7)  and  is  equal  to  the  power 
transferred  across  S by  the  compensating  waves. 


P = Pr  + 


If  equations  (2)  and  (3)  are  evaluated  at  z'  = d+  and 
z 1 - d-,  respectively.  It  turns  out  that  the  powers  are  identical 

T O y 

if  p2(d-)  = n p^(d+).  This  verifies  the  previous  result. 

Figures  18  and  19  contain  curves  of  ■ *o'r~T  as  a function 

2 P2|o>p3 

of  dipole  height,  0z',  for  n = 2.54,  for  wo  thicknesses  of 
dielectric,  0d  = 0.15tt  and  nrr/4.  The  integrals  in  equations  (5-7) 
were  computed  by  numerical  integration.  Comparison  curves  with 
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Pd  = 0 are  given. 


The  discontinuities  at  Pz  = Pd  in  Figs.  18  and  19  are  due  to 
the  above-mentioned  discontinuities  of  the  magnitudes  of  the 
sources*  The  total  powers  are  seen  to  decrease  with  Pz‘  in  the 
dielectric.  The  slight  humps  in  this  region  are  seen  to  occur 
at  the  maxima  in  Figs.  20  and  -21.  For  pz*  > Pd  the  cowers  fluctuate 
around  the  values  for  Pz*  s 0.  Both  converge  upon  ^ as  pz'  ap- 
proaches large  values,  indicating  that  the  coupling  to  the  guided 
modes  and  image  fields  approaches  zero  for  large  Pz'. 

Figures  20  and  21  contain  curves  of  P / P versus  dipole 
height,  pz',  for  the  same  pd's  and  n as  Figs.  18  and  19 o Figures 
20  and  21  are  seen  to  have  maxima  near  the  dielectric  and  conductor 
surfaces  and  to  decrease  uniformly  with  Pz'  and  pz'j^pd.  The 
minimum  between  the  two  maxima  agrees  with  the  minimum  in  Figs.  18 
and  19  for  Pd  = 0.  The  two  maxima  correspond  to  the  two  methods 
of  feeding  a balanced  two-wire  transmission  line:  (a)  a concen- 

trated generator  in  the  terminal  piece,  or  (b)  two  equal  concen- 
trated generators  of  opposite  polarities  in  the  two  wires  at  the 
same  cross  section. 

For  the  smaller  Pz',  the  large  magnitudes  of  the  ratios  of 
Fig.  20  and  21  indicate  that  an  overwhelming  portion  of  the  power 
radiated  by  the  dipole  is  transferred  to  the  guided  or  surface 
waves.  In  these  cases  the  nomenclature  of  compensating  wave  for 
the  radiated  component  of  the  field  is  a preferred  one. 


VIII 

The  Attenuation  Constant  of  a Guided  lode 

Due  to.  the.  Finite  Conductivity  of  the  Ground  Plane 

The  attenuation  constant,  a^,'  for  Pr>>l,  of  the  guided 
mode  due  to  a large  but  finite  conductivity  of  the  ground  plane 
is  derived • 


! 

1 

I 
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FIG.  18  TOTAL  POWER  RADIATED  AND  RADIATION  RESISTANCE  AS  A 
FUNCTION  OF  f3z  FOR  /3d  = 0.  l5Tf  WITH  A COMPARISON  CURVE 
FOR  Bd  = 0 

n =2.54  /?d  = .l5  TT J36  = 0 


FIG.  20  PoUT/P0rUT  AS  A FUNCTION  OF  £z'  FOR  £d  = 0.I5TT 
n2®  2.54  THIS  INCLUDES  CASES  A,B  AND  C AS  &z 
VARIES. 


O 0.2  0.4'  0,6 ySz'/TT  08  10  <-2  1-8 


FIG.  21  ^utAout  ^ A FUNCT,0N  0F  0Z‘  F0R  £d*nTl74, 
n2  s 2.54  THIS  INCLUDES  CASES  A,B  AND  C AS 
£z‘  VARIES 


! 
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Use  is  made  of  the  following  formulas 


1. 


1 atf/a* 


If  the  following  procedure  is  followed  equation  (2)  will 
result  for  cases  A,  B,  and  C. 

a)  Assume  the  fields  are  those  of  equations  (4.1)  - (4.10) 
except  for  a small  component  Erj^(r,o) 

n 

b)  Assume  Pjr  is  that  of  equations  (7.2)  - (7«7). 

c)  Set  H0JG(r,o)  = - lrJ(r,o)  & - Er(r,o)/zg, 

Where?  lrJ(r,o)  is  a quasi-surface  current. 


*s  “ <i-i> 


d)  Compute 


” 2 Ra  / zsl^rJ(r’o)l  rd9 

_ «8 


Hw(r,o)  rd© 


e)  Substitute  e and  b in  equation  1 


3 


2.  a. 


n282  ‘s/’j2-P2t(g2-P~%2)  + n2(\T2  - g2)1 

X.jVSnn*  jp2(n2_i)  + d ( p2-n- V)+n2yV^?) l) 

For  cases  A,  B,  and  C. 

lecture  notes  of  R.  D.  Kodis  and  E.  T.  Kornhauser  at  Harvard 
University. 
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Figure  22  contains  a plot  of  ^2^2  otj  as  3d  varies,  for 
n2  = 2.54,  and.lOn  < 3d  < 1.6rr,  for  the  existing  modes  1 and  2. 

It  is  to  be  noted  that  the  attenuation  constant  for  each  mode 
increases  to  a maximum  and  then  decreases  uniformly  with  in- 
creasing 0d.  These  maximum  values  of  attenuation  appear  approxi- 
mately halfway  between  the  cutoff  thicknesses.  The  attenuation 
constants  for  the  higher  modes  are  lower  than  those  for  the  lower 
modes  since  the  factor  d appears  in  the  denominator  of  (2)  and 
the  higher  modes  appear  with  larger  3d. 


IX 


The  Radiation  Resistance 

The  radiation  resistance  may  be  computed  by  referring  it 
to  the  uniform  current  on  the  dipole,  I(z'). 

2P 


(1) 


R 


e _ 


|l(z*)l; 


where  I(z')  = - for  both  the  Abraham  and  Hertz  dipoles 


and 


Hence: 


2h 

h = the  half-length  of  the  Hertz  dipole  = the  length  of 
the  Abraham  dipole. 


(2) 


e 


R ” = 4tteP 

240 32h2  a>33  |p(z‘)  ft 


This  is  plotted  with  P in  Figs.  18  and  19 • 
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Appendix  A 

To  illustrate  the  type  of  integration  performed  to  yield 
equations  (3*5),  (3*6),  (3.7)  and  (3,9),  consider  the  integral 
in  (2.29),  after  it  has  been  transformed  by  (3.1),  to  be  inte- 
grated over  of  Fig.  3 • 


, . 4(z+z‘-2d)  , 

(A-l)  _ I - i I ^ 5T H/CKrJMA. 

sC+mn  canh  md 

It  is  convenient  to  make  the  following  substitutions 
(A-2)  R',sin91  = z + z*  -2d 5 R^osO-j^  = r Fig.  2a 

(A-3)  ^ = Psina  5 yi2  - \2  = pcosa 


The  Riemann  sheet  of  interest  in  the  \-plane  of  Fig.  3 
transforms  into  the  curved  strip  in  the  a— plane.  Fig.  A-l. 

The  transformation  of  the  branch  cuts  (which  are  not  branch  cuts 
in  the  o-plane)  are  defined  by  sin  x cosh  y=T...  Under  the 
assumption  that  0R“  cosfij^l  so  that  H0*  may  be  replaced  by  its 
asymptotic  form  for  large  argument,’  and  substituting  (A-2)  and 
(A-3)  in  (A-l),  (A-l)  becomes s 


I * - 
(A-4) 


« 


-in/4e2 

4n 


M 


2s  Ina 


nR" {Jcoa&L 


"ipRw(  cosasinO-j+sinpcosO-,) 

..  gQg.ftQ - . . 

[ i0cosa+0n”Vn2-s  in*fc  tan(  0d  Vn2-sin&  ) ) 


This  is  of  the  proper  form  for  an  integration  by  the  method 
of  steepest  descents,*0*  *2  when  R">>d* 

The  saddle-point  is  at  a " n/2 -9. 

The  path  of  steepest  descent,  W,  is  defined  by 
sin(x  + ©2^cosl1  y = 1 (Fig’  A-l). 

If  the  assumption  is  made  that  the  poles  are  not  near  enough 
to  the  a = tr/2  to  affect  the  integrand  before  the  exponential 


* 

•t 

t 
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factor  has  attenuated  it  to  negligible  value  ’ then  the 
integration  over  W results  ini 


iPR"  sin©., 

( A- 5)  - i — ....  

ipsin©1  + pn~^  /^ln2- ccs2©^tan(  Pd'^-cos2©^ 


and  it  may  be  observed  that? 

da  + 2nl^>  ' Residues  enclosed  by 


(A-6)  I = I . . • do  “ " f*  • • » 

J jr„ 

W,  W 


W and  W^. 


When  W passes  through  a pole,  one-half  the  residue  is 
included^’  as  a contribution  to  the  integration*  If  the  pole 
is  included  between  W and  W^,  the  whole  residue  is  contributed 
to  (A-6),  while  if  the  pole  is  excluded  by  path  W-W,9  there  is 
no  contribution  from  it* 

It  Is  relatively  easy  to  verify  that  W passes  through,  a 
pole  when  ©1  = cos^p/kj,  and  from  this  relationship  the 
coefficients  Aj(O^)  °**  the  residue  waves  arise.  The  residue  terms 
are  evaluated  by  well-known  methods. 

The  rapid  decrease  of  the  coefficient  Aj(©^)  at  Its  critical 
value  is  due  to  the  fact  that  for  large  PR  the  path  passes  very 
rapidly  through  the  pole  as  ©^varies. 

The  results  obtained  thus  far  have  been  restricted  to  PR" 
cos©^>>l  by  virtue  of  the  assumption  that  the  Hankel  function 
could  be  replaced  by  its  asymptotic  form.  Actually,  the  re- 
striction is  less  strict  for  a good  approximation  and  it  may  be 
modified  to  PR**»1  and  (PR"cos©1)  >>1.  This  leaves  only  a 
small  region  near  ©^  = n/2  where  the  results  developed  are  not 
valid.  The  restrictions  may  be  further  modified  by  the  followings 


See  Appendix  B for  the  exact  condition. 


I 

| 


s 


FIG.  A-l  PATH  OF  INTEGRATION  IN  THE  oL=sin"'  x/3  PLANE 
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Consider  (2.29)  for  r = o,  0^  = tt/2,  3R“  = 3z>>! 


(A-7) 


^-'tC z+z • -2d ) 

i •« + 


_p 

mn  tanh  md 


XdX 


Performing  the  same  transformations,  (A-2)  and  (A-3),  and 
integrating  over  half  the  path  W (the  part  on  and  below  the 
x-axis  of  Pig.  A-l)  in  the  emplane,  the  following  is  obtained: 


(A-8) 


i e ip(z+z»-2d) 

2™  (z+z '-2d)( i+n”  tanpnd), 


where  there  are  no  residue  terms  since  3z  >>1  and  ©^ 
cos^P/XjCCAj^©)  = 0). 


It  is  seen  that 


I 


r*=o 


lim  It 

®1=  n/2  | eq.  A-6 


n/2  > 


It  is  therefore  a reasonable  assumption  that-equation  (A-6) 
is  valid  for  all  ©x,  0<91  <tr/2,  0R»>>1. 


The  integrations  of  equations  (2-30),  (2.31),  and  (2.33)  are 
performed  in  a similar  manner. 


jommMUfciiittitrti — - — »■  wti  -iMiin—— 
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Appendix  B 

As  an  illustrative  example  of  the  application  of  the  method 
of  B.  L.  Van  Der  Waerden^  to  the  integration  of  equations  (2.32) 
and  (2.34),  consider  the  following  integral,  obtained  by  applyiing 
(3.1)  to  (2.32). 


(B-l) 


00 

-4; 

-CD 


(4-mn“2)coshmz 'coshmz  e “in<5?-.  B1 


o (Xr)dX 


m (■$  coshmd  + mn“2  sinh  md) 


Since  the  integrand  of  (B-l)  has  two  extra  branch  cuts, 

(3)  and  (4)  in  Pig.  3,  at  X = + n0,  and  one  of  these  intersects 
the  path  of  steepest  descent  for  the  integrand  in  (B-l)  if  the 
method  of  steepest  descents.  Appendix  A,  is  used,  the  integration 
by  the  method  of  steepest  desents  is  complicated.  The  Van  Der 
Waerden  method  contains  a systematic  procedure  for  integrating 
around  branch  cuts. 

In  order  to  put  (B-l)  in  Van  der  Waerden* s form  the  following 
transformation  is  performed. 

(B-2)  X = in 

Substitute  the  asymptotic  representation  for  large  Xr  for 
H^CXr). 

The  transformed  \i  = - IX  - plane  is  shown  in  Pig.  B-l  where 
(1),  (2),  (3),  (4),  and  (5)  are  the  branch  cuts  that  appear  in 
the  X - plane,  Pig.  3* 

Then* 

(B-3)  I * 

Where  r corresponds  to  Van  Der  Waerden* s X* 

In  the  n-plane  the  saddle-point  of  the  ^-plane  becomes 
another  branch  point. ^ In  this  case  there  is  no  added  complication 


/ 


W e"r^dv 


for  0r  >>  0d, 


\ 


I 


v «*&•**;  tt&i  "**»  ■•  * •—■***» 
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because  this  added  branch  point  is  coincident  with  the  one  at 
M-  =-ip. 

For  large  r and  for  a contour  shifted  far  to  the  right  in 
Fig.  B-l,  the  integral  reduces  to  integrations  over  the  branch 
cuts  and  poles.  The  former  yield  rapidly  converging  series  for 
large  Xr  and  for  poles  not  to  close  to  the  branch  points. 

Performing  the  integrations  exactly  as  Van  Der  Waerden  does, 
the  following  is  obtained: 

k 7 

(B-4)  I = 2fti  ^ Residue  - ^ e in^r+  0(r”2),' 

J=1 

—2  —2 

where  0(r  ) indicates  a term  of  order  r and  yields  an  estimate 

of  the  error.  


The  residue  are  computed  in  the  usual  way  and  are  modified 
by  equation  (3.2)  to  yield  the  summation  terms  of  (3.8). ~ 

If  in  equation  (2.32)  it  is  assumed  that 


0z* 

and 

Pz 


< Pd  «0r, 


then  0RV  * PR'  = Pr  and  equation  (3*8)  is  obtained. 


It  has  been  assumed  that  the  poles  are  not  too  near  to  the 
branch  points.  An  estimate  of  not  too  near  is  a result  of  the 
Van  Der  Waerden  method.  It  Is: 

-rlS-fcjl 


(B-5) 


a) 

b) 


and 

-rfnp-xj 


«'0(r"2)«0(r"1) 


If  the  integrations  in  Appendix  A had  been  performed  by  the 
method  in  Appendix  B,  the  definition  of  the  pole  being  not  too 
near  the  branch  point  would  be  found  to  be  (B-5a).  Hence  equation 
(B-5a)  and  (B-5b)  are  necessary  restrictions  for  all  the  equations 

(3.5)  - (3.10). 
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